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Abstract 

Following a strategy introduced earlier by the authors, we show that it is 
possible to extract \V u b\ from the cut hadronic mass spectrum of B decays 
without large systematic errors which usually arise from having to model the 
Fermi motion of the heavy quark. We present a closed form expression for 
|Fub|/|Vt s | which is accurate up to corrections of order a 2 s , a s p, (A/mj,) 2 / p, 
where p is the experimental cut (s c /m 2 ) on the hadronic mass used to veto 
charmed decays. Modulo duality violation errors, which are intrinsic to all 
inclusive predictions, we estimate the theoretical error in the extraction to be 
at the 5% level. 
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I. INTRODUCTION 



An accurate determination of V u b would fill a gaping hole in our understanding of the 
Cabibbo-Kobayashi-Maskawa (CKM) sector of the standard model. Indeed, while it is 
not clear that we will be able to cleanly determine the angles of the unitarity triangle, 
an accurate measurement of \V u b\ could determine whether or not the CKM description is 
correct by measuring the lengths of the sides. 

While in principle one could measure \V u b\ quite easily in a systematic approximation 
to QCD [T|] from an inclusive measurement of the uncharmed semi-leptonic branching ratio, 
in practice, it seems to be much more difficult. The problem lies in rejecting the charmed 
decay background. This cut causes the canonical expansions in a s and A/m& to breakdown. 
The imposition of the cut changes the expansion parameters to a s log p and A/ (mj, p), where 
p parameterizes the cut and is numerically small. It may be the case that one or both of 
these parameters can be of order 1, thus necessitating a reorganization of the expansion. 
Here we will not review the systematics of these expansions, as they have been thoroughly 
discussed in the literature || . 

The nature of the abovementioned series depends on the kinematic variable with which 
we choose to cut. Each choice has advantages and disadvantages. The cut lepton energy 
spectrum, E\ > (M| — Mq)/(2Mb), is relatively simple experimentally since there is no 
need for neutrino reconstruction. Theoretically, this choice of cuts is challenging due to the 
breakdown of the non-perturbative series in A/(m& — 2Ej). Furthermore, the fraction of the 
b — > u transitions included after the cut is only about 10%. The cut hadronic invariant mass 
spectrum is more challenging from an experimental viewpoint, but has its virtues ||. First 
of all, while this observable smears over the same range of hadronic masses as the cut lepton 
spectrum, it is weighted more towards states of larger invariant mass. We thus would expect 
local duality to be more effective in this case. Furthermore, this cut rate includes around 
40 — 80% of the b — > u transitions ||. In both of these cases the non-perturbative series 
must be resummed. The necessity for resummation of the perturbative series can only be 
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determined a posteriori. For example, in the case of radiative decay with an experimental cut 
of 2.1 GeV on the photon energy, it was found || that the resummation of the perturbative 
series was not necessary. In |7j it was shown that, ignoring the non-perturbative structure 
function, the perturbative resummation is crucial, and that the cut hadronic rate can not 
be calculated in a systematic way. However, as was emphasized in |7|], this conclusion may 
not stop us from extracting \V u b\ since the inclusion of the structure function can greatly 
soften the effects of perturbative resummation. 

In [U it was shown how the endpoint data from radiative decay can be used to eliminate 
our ignorance of the Fermi motion of the heavy quark (i.e., the non-perturbative structure 
function) in the electron end point spectrum. In this paper we will apply these ideas to 
the hadronic mass spectrum. The next section is a synopsis of the results given in for 
the resummed hadronic invariant mass spectrum, which is followed by a section dedicated 
to applying the formalism of 0] to derive a closed form expression for \V u b\ in terms of the 
hadronic mass spectrum. The final section of the paper includes the results and a discussion 
of the errors. 



II. INFRARED FACTORIZATION AND RESUMMATION 

The factorization and resummation for inclusive B decays in leptonic variables are dis- 
cussed in Ref. ||. It was extended to hadronic variables in Ref. [[/]]. Here we will briefly 
review the results. Consider the inclusive semi-leptonic decay of the b quark into a lepton 
pair with momenta q = (pi +p v ) and a hadronic jet of momenta pt- We define the following 
leptonic kinematic variables 

Vo = , 1 

m b 

X = (3) 
m b 

where v = (1,0), and partonic kinematic variables 



h = (5) 

In terms of the leptonic variables, Sq = (1—yo+y) and h = 2 — y , one can see that in the end 
point region of the electron energy spectrum when x — > 1 with y < 1, the invariant mass of 
the jet approaches zero with its energy held fixed. In addition, the jet hadronizes at a much 
later time in the rest frame of the B meson, due to the time dilation. Factorization exploits 
this and separates the particular differential rate under consideration into subprocesses with 
disparate scales. This factorization fails when the jet energy vanishes in the dangerous region 
y — > x — > 1. However, this problematic region of phase space is suppressed because the rate 
to produce soft massless fermions vanishes at tree level. 

In terms of yo and y, the triply differential rate, which factorizes into hard, jet and soft 
subprocesses f8j, may be written as 

1 d 3 T 



T dy dy dx 



dzS(z) m b J[m b h(z - £)> A 4 ] H(m b h/fi), (6) 

where x m = Ms/nib, n is a factorization scale, and £ = (l — y)/(2 — y ) is analogous to the 
Bjorken scaling variable in deep inelastic scattering, z — 1 + k+frrit,, where k + is the heavy 
quark light cone residual momentum. S(z) essentially describes the probability for the b 
quark to carry light cone momentum fraction z and allows for a leakage past the partonic 
endpoint, as can be seen explicitly in the upper limit of z. 

It is convenient to change variables from y to £ and perform the x integration to yield 

^^7 = 2(2-t/ ) 3 (2yo-l) dzS(z)m 2 b J[m 2 b (2-y )(z~O}H[m b (2-y )}. (8) 
i o ay at; J$ 

To proceed, we take the iVth moment with respect to £ in the large iV limit. In the region 
s ~ and z ~ £ ~ 1, one can replace J\m L h{z — £)] in Eq. (g) with J[m 2 h{l — £/z)]. This 
replacement is permissible to the order we are working. We then obtain 



M N = d£ £ 
Jo 



,v_i 1 d 2 Y 
T df dy 



= 2(2 - 2/0 ) 2 (2y - 1) S N J N [m 2 b (2 - y )//i 2 ] H[m b {2 - y ) /fi] + 0(1/N), (9) 
Mml/fi 2 ) =m\ f'dyy^JimKl-y),^, (10) 



5^= / dzz N S(z). (11) 
•/o 

The soft moment Sn further decomposes into a perturbative piece, which accounts for soft 
gluon radiation and a non-perturbative piece which incorporates bound state dynamics and 
serves as the boundary condition for the renormalization group equation M. 

S( g )= r m ^f[ mb ( s -l)]a(z/s), (12) 

Jz S 

Sn = In cat, (13) 

where f(s) = (B(v)\ b v 5(s—iD + ) b v \B(v)) is the non-perturbative structure function defined 
in Ref. [f|, which resums an infinite number of operators all in the same order of A/m&. A 
similar expression holds for the inclusive B — > Xg'j in the end point region of the photon 
spectrum || 



Tq Jo 77 dx y 
= S N J N H\ (14) 
G 2 F \V*V tb \*aCtml 
10 = 32^ • (15) 



Here C-j is the Wilson coefficient of the 0-? operator evaluated at the scale m b |L0| and 
x 7 = 2E 1 /m b . 



In moment space the soft and jet functions [|Tl|] have been calculated to next-to-leading 
(NLL) logarithmic order []2| and are given by 

a N J N = exp[log(iV)#i(x) + gj(x) +9si(x,yo)) (16) 
a N Jl = exp[log(A0<7i(x) + gl(x)l (17) 

where \ — a s( m l)Po l°g N, and gi, g\ and g s \ can be found in [Q. The hard parts are given 
in R: 



H(yo) = l 2n 



3tt 



, o 8- 10y ,__, n \ ( nT * / _ . , 5 , 2^ 



4 log 2 (2 - y ) + -= ^log(2 - y ) + 2U 2 (y - 1) + - + 



2J/0-1 2 3 



(18) 



To get back the physical spectra from the moment space, the inverse Mellin transform 
has to be evaluated at NLL accuracy as well. If we are willing to ignore the effects of the 
structure function we can get a closed form expression for the resummed rate. To this end, 
we use the results derived in the Appendix of Ref. ||] 
1 d 2 T 



T dy d£ 



2(2-y o y(2y ~l)H(y ) 

d f gZgi(a s A)0+92(a! s /3oO ~| 



where I = — log (—log £) ~ — log(l — £) and g 2 = gj + g s i- The ^-function defines the 
differential rates in a distribution sense, as rj — > 0, and turn the singular terms into the '+' 
distributions. 

This result would only be useful in making a physical prediction if we had a measured 
structure function with which we could perform a convolution. Instead, following J7| we will 
take the inverse Mellin transform of the ratio between the moments of the semi-leptonic and 
radiative decay rates. 



III. EXTRACTION OF \V ub \ 



The key observation to extract \V u b\ is that the soft function, Eq. (|12]), is universal in 
inclusive B decays. Using ( |Tjj ) - ( p]) we can rewrite (^) as 



M N = 2(2-y y(2y -l) 



I) E^A e 9si(x,yo) M 



(21) 



Taking the inverse Mellin transform we arrive at 
Vl £T : = 2(2 - y )\2y - 1) 

dT^ ( 



T dy d£ 



x / H(y )du- 
J$ du 



dz 



9(1 — z — rj)e 



(22) 



where 

H{y ) 



H(y ) 



Hi 
1 _ 2a * 



3tt 



4 log 2 (2 - y ) + 8 ^ log(2 - y ) + 2U 2 (y - 1) - 4 



(23) 



4 

0d(x, 2/o) = 5-5- log(2 - y ) log(l - x), (24) 

% 2 = — a s fio log(l — 2) and 2 = We have now rescaled to "physical" variables which 
are normed to the meson mass Mb- We have also dropped a factor of A/rrib which would 
multiply yo in these renormed variables, as it would only be relevant in the region where 
factorization breaks down, yo ~ 2, which, as previously mentioned is the suppressed part of 
phase space. 

Now to get the cut integrated rate we change the order of integration so that integral 
over photon energy is the last to be performed. This leads to a result which can most easily 
be integrated over the data for the radiative decay. We find the following expression for the 
rate with a cut c = s c /M B on the hadronic invariant mass spectrum: 

(fTf~ f 1 Muz— 2 rl /•-rrp-J Muz— 2 

du /wcZj dz J 2 _^_ dy ° + 1^ du U-Jn dz J 2 _^_ d V° 
2 u 1 — uz 1 + c u 1 — uz 

+ I du I dz i dyoj-pf I(u, y Q , z) (25) 



l+c u(l+e) 



I(u,y ,z) = 2(2 - yo )\2y -l)H(y )^- (-z-f) \9(1 - z - V ) e^^] . (26) 



du y dz 

We can estimate the effects of resummation on the cut rate by using a model for the 
structure function. In Figure 1, we show the effects of the resummation of the end point 
logs as a function of the cut using the same model as in ||. We see that the effects of 
resummation are small, and thus the end point logs do not form a dominant sub-series, and 
the unknown piece of the two loop result dominates the perturbative uncertainties. Note 
that this is not inconsistent with the results of where it was found that the effects of 
resummation were large, as the structure function mollifies the effects of the resummation. 

If we now expand g s i to order a s , we can write down a simple expression for \V u b\ as a 
convolution of the data from radiative decay. We may use the approximation 
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FIG. 1. Percentage change due to resummation. The solid line is the fully resummed result. 
The dotted line is the expansion of the resummed result to order a s , while the dot-dashed line is 
the expansion of the resummed result up to order . 



^(6(l-z)e 9 ° l )n5(l-z) + 
dz 



Aa s ( log(2 - y ) ' 



3tt 



1 - z 



We then find 



ar( c ) = -4[i (c) + / + (c)], 

1 n 



i 

l+e 



4u-2 



dT"< ~ 



/o(c) = du L— dyo + i^- du JL* dyo i 2(2 ~ 2/0)2(22/0 _ 1 ^ H( - y °^ 

\ 1 l — u 1 + c u 

J+(c) = JduJdzJ dy 2(2 - y ) 2 (2y - l)z 



4u-2 

%oj2(2-y ) 2 (22/o-l)-^ 

u / 

dT^Aa s Aog(2-y ) 



37T 



1 - z 



The plus distribution is defined as 



1 — x Jo 1 — a; 



(27) 



(28) 

(29) 
(30) 



(31) 



The regions of integration in Eq. ( |3"0"D are the same as those in Eq. (P5|). The integrals 
over y and z can all be done analytically. The final result is rather large and complicated, so 
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we choose not to include it in print. The resulting expression is an integral which involves 
a known function of u and c and the radiative decay endpoint data. Thus we may write 

\V uh \ _ \ 6aC 7 (m b )HT{c) Y 

\V ts \ \ n[I Q (c) + I + (c)] J ■ 1 ] 

The dominant source of error, modulo the usual uncertainties inherent in all inclusive pre- 
dictions, will come from unknown effects at order or s , and A 2 /(cmj), each of which should 
contribute only at the few percent level. 



IV. RESULTS, CONCLUSIONS AND OUTLOOK 

The extraction of \ V u b\ using the cut hadronic mass spectrum has the theoretical advan- 
tage over the lepton spectrum since the mass cut includes a larger fraction of the b — > u 
transition. In addition, we expect local parton-hadron duality to work better since more 
resonances will contribute to the cut mass spectrum. The real issue for the viability of 
this method is the resolution. In particular, how good will the "best" resolution be for the 
invariant mass? This will be limited by the inability to detect neutrals, as well as the charge 
particle detection inefficiencies. Presently, the resolution is not good enough to eliminate the 
charmed background, which has to be modeled. Hopefully, it will be possible to eventually 



push the cut further |[L2| . If the resolution does become good enough to really eliminate 
the charmed background, then there is an additional option available. It has been pointed 
out that the leptonic mass spectrum is effectively less sensitive to the Fermi motion of the 



b quark ]T3[. Recently, this fact has been used to determine an expression for \V u b\ in terms 



of the cut lepton mass spectrum, without needing to use the data from radiative decay JTJ 



However, the resulting expansion parameter is A/Md, and grows as the cut is raised above 
(Mg — M^) 2 . The authors in [TJ] calculated the leading non-perturbative correction and 



found it to be anomalously small. They made an estimate for the sub-leading corrections 
and found that they grow rapidly as the cut is raised from its smallest possible value. 

In the end it is clear that no one single extraction of \V u b\ should be trusted, given 
the fact that our expansion parameters are never as small as we would like. We will only 
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gain confidence in the extractions after there is convergence among several independent 
extractions. 
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